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The transverse spin structure of the pion at short distances 
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Abstract: We study the form factors of the quark tensor currents in the pion at large squared 
momentum transfer Q^. It turns out that certain form factors can be evaluated using coUinear 
factorization, whereas others receive important contributions from the end-point regions of the 
longitudinal quark momenta in the pion. We derive simple analytic expressions for the dominant 
terms at high and illustrate them numerically. 



1 Introduction 

The structure of the pion at short distances unites two 
characteristic features of quantum chromodynamics. 
On the one hand, the pion plays a unique role among 
hadrons as the Goldstone boson of spontaneous chiral 
symmetry breaking. On the other hand, asymptotic 
freedom is central for understanding its structure at 
short distances, where quarks and gluons interact per- 
turbatively as in any other hadron. Moreover, many 
studies of hadron structure are very much simplified 
when one deals with spin-zero hadrons, and the pion 
is probably the spin-zero hadron for which most quan- 
titative information is available, both from experiment 
and from calculations in lattice QCD. A versatile tool 
to describe hadronic structure is given by generalized 
parton distributions or, equivalently, by the form fac- 
tors of a tower of local quark-gluon operators contain- 
ing an increasing number of covariant derivatives. 

A perhaps surprising feature of the pion is that is 
has a non-trivial spin structure. An instructive quan- 
tity to describe this structure is the distribution p{x, b) 
of quarks with longitudinal momentum fraction x and 
transverse distance b from the center of the pion [1 . 
Due to parity invariance this distribution cannot de- 
pend on the longitudinal quark polarization. However, 
the distribution of quarks with transverse spin s has 
a polarization dependent part, which is proportional 
to (s X 6)^ and was found to be sizeable in a recent 
lattice study [2]. This polarization dependence can 



be quantified by the form factors of the quark tensor 
operator qia"^ q and its analogs containing covariant 
derivatives. The present work is concerned with these 
tensor form factors at high momentum transfer, or in 
other words with the correlation between the trans- 
verse polarization and the transverse position of quarks 
very close to the center of the pion. 

Form factors at high momentum transfer have 
played a key role in the early development of meth- 
ods for calculating exclusive observables in QCD [3l|4]. 
They continue to provide an important area for apply- 
ing factorization, with close links to the physics of ex- 
clusive B meson decays. In the limit of infinite momen- 
tum transfer form factors can be described within 
standard collinear factorization, but extensive studies 
of the electromagnetic pion form factor Ft^{Q'^) indi- 
cate that at experimentally accessible values of this 
description receives important corrections, see for in- 
stance O [6l [ll EJ [9] . In the present work we aim at pro- 
viding a baseline for the large behavior of the pion 
tensor form factors i?Tnj , and we will use a very simpli- 
fied extension of the collinear factorization framework 
that allows us to obtain expressions in compact analyt- 
ical form. We do therefore not expect our results to be 
quantitatively reliable at moderately large , and we 
will in particular refrain from comparing to the lattice 
calculations in [5], which go up to « 2.5 GeV^. On 
the other hand, our analytic expressions may be of use 
if one wants to devise parameterizations of BTndQ^) 
that have the correct behavior at large . 
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The large behavior of pion tensor form factors 
is also interesting because it involves pion distribu- 
tion amplitudes of twist three, which have a particu- 
lar behavior at the end-points of the momentum frac- 
tion variable [TU]. We find that for certain form fac- 
tors Bxni the formulae obtained by using collinear fac- 
torization have end-point divergences and hence need 
to be modified. This is similar to other cases where 
twist-three pion distribution amplitudes appear, such 
as spectator interactions in exclusive B decays [TTl [13] , 
pion electroproduction ep — > e7r+n with transverse po- 
larization of the exchanged virtual photon [14j . and 
certain power corrections to Ft^{Q^) [T51 116) . 

This paper is organized as follows. In the next sec- 
tion we set up the calculational framework used in the 
present work. In Sect. |3] we extract the contributions 
from the hard-scattering graphs that dominate in the 
large limit and derive simple analytic expressions 
for the form factors Btui- In Sect. 2] we present some 
numerical illustrations of our results, and in Sect.[5]we 
summarize our findings. 

2 Setting up the calculation 

The tensor form factors of the pion parameterize the 
matrix elements of the local operators 



T A S qia 

{a,Pi) (/Ji,...,/3„) 



(1) 



where 3^ = d^-igA^ with cf'^ = ^0/^ -^1^) is the co- 
variant derivative. Here S and A respectively denote 
symmetrization and antisymmetrization in the indi- 
cated indices, and T denotes the subtraction of traces 
in all index pairs. These operations, which project 
on operators with twist two, can be implemented in a 
simple way by contraction with two constant auxiliary 
vectors a, b satisfying = ab = and &^ 7^ [17]. 
The tensor form factors arc then given b}0 

(7r+(p') I uia^'^ac^bp {iDa)''-^] ti+ {p)) = (aP)"-'^ 



i = 



with = —{p — pY and 



P=^{p + p'), 



e-^^- (3) 
a{p + p') 



The form factors in ^ refer to u-quarks; those for 
c?-quarks follow from isospin symmetry and read 



B^,„ = (-l)"i?. 



Tni- 



(4) 



The form factors can be written as Mellin moments of 
generalized parton distributions of the pion as shown 
in |17| . but we will not need this representation here. 

In the collinear factorization formalism and at lead- 
ing order in the matrix element ^ receives con- 
tributions from the graphs in figure [T] Due to the 
covariant derivatives, the operator ([l} contains terms 
with zero to n — 1 gluon fields. Graphs (a) and (b) 
correspond to the term without gluon fields, i.e. to 

uia''^aabf3{iday'-^u (5) 

in ((2). The same graphs describe the electromagnetic 
pion form factor if one inserts the electromagnetic cur- 
rent instead of the current in ([5]). Graph (c) corre- 
sponds to the terms in ([Ij that have exactly one gluon 
field, i.e. to 



n-l 

E 



uia"^ Gabp (ida) 



n-l-j 



{gAa){iday'^u (6) 



in Terms with more than one gluon field do not 
contribute at this level. 

When calculating the hard-scattering part of the 
graphs we neglect the pion mass, so that the pion mo- 
menta p and p' are purely lightlike. We use them to 
define the two light-cone directions required for speci- 
fying the distribution amplitudes of the pions, working 
in a reference frame where the incoming pion moves in 
the positive and the outgoing pion in the negative z 
direction. As indicated in the figure, we write the u- 
quark momentum as up + k in the incoming Tr~^ and 
as vp' + k' for the outgoing 7r+, with the light-cone 
momentum fractions u and v ranging from to 1. 
The vectors k and k' are transverse to both p and 
p' . We neglect the small momentum components of 
the quarks and antiquarks, i.e. the component along 
p' in the incoming pion and the component along p 
in the outgoing one. Note that {up + k)"^ = k"^ and 
{vp' + fc')^ = fc'^ are in general not zero — we will com- 
ment on this shortly. 

Since the tensor operators ([T]) have odd chirality, we 
need one chiral-even and one chiral-odd pion distribu- 
tion amplitude in the graphs to obtain a nonvanishing 
hard-scattering amplitude. Since there is no chiral-odd 
pion distribution amplitude with twist two, we must 
go to twist-three level. The relevant distribution am- 
plitudes have been introduced in [10] . After a Fourier 
transform from the position representation used in [10] 
to momentum space, the projection operators for the 
incoming and the outgoing pion respectively read [11] 



A factor i is missing on the r.h.s. of eq. (71) in [Tf] . 



i)a{u) Pap'i3 ^ ^^^^ d 

du pp' dk^ 



2 



p -p 



up + k I 
-®- 



vp' + k' 



(a) 



up — k 



vp' - k' 



P -P 



up + k 



4 



I vp' + k' 



(b) 



up — k 



vp' — k' 



p -p 



up + k I 



vp' + k' 



4 



(c) 



up — k 



vp' - k' 



Figure 1: Graphs for the matrix element ([2]) in the 
Hmit of large Q^. The crossed circle represents the 
insertion of the relevant current operator, given by ([5]) 
for graphs (a) and (b) and by ([6]) for graph (c). The 
blobs stand for the sum of twist-two and twist-three 
distribution amplitudes as specified in ([7]). 
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In ([U the pion mass can of course not be neglected 
since one is dealing with a non-perturbative quantity. 
For the twist-three distribution amplitudes we take the 



asymptotic forms under evolution, 

(t>p{u) = 1, 4>a{u) (9) 

where here and in the following we use the notation 
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(10) 



The normalization constant f^-^ associated with the 
twist-three quark-gluon-quark distribution amplitudes 
of the pion asymptotically evolves to zero TOl. In the 
limit where <j)p and <j)„ take the form ([5]), the graphs 
in figure [T] therefore give the full answer for the matrix 
element ^ . Conversely, the consideration of distribu- 
tion amplitudes deviating from ([SJ would require the 
inclusion of graphs with an additional gluon in one of 
the pion distribution amplitudes and thus considerably 
complicate the analysis. Since in this work we aim at 
understanding the basic behavior of the form factors 
at large Q^, we consider the restriction to the asymp- 
totic forms (j9l) to be sufficient. On the other hand, we 
can easily keep the general form 

4>{u) = 6uug{u) (11) 

of the twist-two distribution amplitude, where 



(12) 



is the usual expansion in Gegenbauer polynomials, 
with coefficients a„ that evolve with a simple multi- 
plicative factor at leading order [3l [J. With ([9]) to 
the factorization scale dependence of the projec- 
tors ([7]) is then given by 



4//3o 



, , , ,faMY''-^''° ,,,, 
a„(^) = a„(/io) [j^~(J^j' J (13) 

at leading logarithmic accuracy, where as(/i) is the 
one-loop running coupling, /3o = 11 — 2ni?/3, and the 
first few anomalous dimensions read 72 — 50/9, 74 = 
364/45, etc. The scale dependence of /i^ simply reflects 
the running of the quark masses in ([5]) . 

An alternative form of the projector ([7]) was derived 
in section 3.2 of [H], which had earlier been used in 
[151 116) . This derivation requires one to keep the small 
components of the quark and antiquark momenta in 
the intermediate stages of the calculation and to ad- 
just them such that both the quark and the antiquark 
attached to the pion wave function are exactly on shell. 
Having the external quarks and antiquarks of the hard- 
scattering subprocess exactly on shell is certainly an 
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attractive feature of the calculation, especially from 
the point of view of gauge invariance. It comes, how- 
ever, at the price of violating momentum conservation. 
Consider for definiteness the quark and antiquark mo- 
menta in the incoming pion: 



Up + k + Wqp', kq 



up 



WqP 



(14) 



For generic values of u and k one cannot have both 
fc^ = fc? = and Wq + Wq — Q (for this it does not 
matter whether one neglects the pion mass or not). 
In our calculation, we choose to be consistent with 
momentum conservation neglect the small components 
Wqp' and Wqp' . We will explicitly check that gauge 
invariance holds for the class of covariant gauges and 
within the accuracy of our calculation. 

As explained in [TT] , the derivatives with respect to k 
and k' in the projector ([7]) act on the hard-scattering 
kernel before one takes the coUinear limit by setting 
k = k' = Q. However, we will see that for some of 
the form factors -By„j the collinear limit cannot be 
taken since the integrals over u and v diverge at their 
end-points for fc = fc' = 0. To keep the intermedi- 
ate steps of our calculation well-defined, we introduce 
transverse-momentum dependent factors S(w, fc^) and 
I](w, fc'^) for the incoming and outgoing pion. These 
factors are real- valued and normalized as 



(fk T.{u,k^) = 1. 



(15) 



In a more sophisticated approach, which has for in- 
stance been used in jl4j , one would multiply the differ- 
ent terms in $ and with different factors and inter- 
pret the result as pion light-cone wave functions that 
depend on both a longitudinal momentum fraction u 
or V and on the transverse parton momentum. Fur- 
thermore, in the spirit of the modified hard-scattering 
approach, one should include Sudakov factors for each 
pion, which resum a class of large logarithms from 
higher-order corrections and depends on the momen- 
tum fractions, the transverse parton momenta and the 
hard scale Q in a non-trivial way 16,. Formally, the 
Sudakov factors alone would already remove the end- 
point divergences of the u and v integrals, but for a 
wide range of hard scales the resulting integrals 
will receive large contributions from phase space re- 
gions where parton virtualities are low and the pertur- 
bative expression of the Sudakov factors is not justi- 
fied (see [18] for a detailed analysis of the situation in 
semileptonic B tt decays). Moreover, even a calcu- 
lation with Sudakov factors but without a nonpertur- 
bative transverse-momentum dependence of the pion 
wave function would not readily yield simple analytic 
expressions. Since the latter is what we are aiming 
for in the present work, we will use a global factor 



S(m, fc^) S(w, fc'^) as a minimal version to regulate the 
intermediate steps of our calculation and simplify the 
resulting integrals in the end, see eq. ([55)1 below. 

With these preliminaries we can write the large- 
limit of the matrix element we are interested in as 

{tt+{p') I uia°''^aabi3 {iDa)"-^u \ TT+{p)) 

— 4:TTas^^ 6f^ I^T^ J dudv dPkSk' Yjiu^k^) 

X E(u,fc'2)/(u,w;fc,fc') (16) 

with 
f{u,v;k,k') 



1 



X Tr$ 



D 



XfJ, 



uvQ"^ - (fc - k'Y 



I / \n— 1 -^ap - 

+ [aty) la ^ ■ 



■r 



n-l 



Gab/S , (17) 



where the last three lines of ([T7| respectively corre- 
spond to graphs (a), (b) and (c) in figure [T] The fac- 
tors 

alu — ^{u ~ u)ap + ^ap' + ak = aP{u — ^u) -|- afc, 
a^-u — 5 (^^ ~ v)ap' + ^ap + ak' = aP{v + ^v) + ak' 

(18) 

come from the derivatives i^a = ^ {id — id)a in the 
operators (O and The denominator of the gluon 
propagator in all three graphs is uvQ'^ ^ {k — k')"^, and 
the quark propagators in graphs (a) and (b) have de- 
nominators uQ^ — k^ and vQ^ — fc'^ . Note that we are 
using a Minkowskian scalar product for the vectors k 
and k' , so that /c^, k'^ and (fc — fc')^ are negative. In 
Feynman gauge, the numerator of the gluon propaga- 
tor is Dx^ = g\fj_ and the fermion trace evaluates to 



f{u,v;k,k') 

X {aLr-' 

a 



fl+f2 + ^{f" + f"-f") 
1 



U 



uQ^ - fc2 toQ2 _ (fc - k'y 

n n-l 



1 



^(a;„)^'-i(a/„)"-i 



1 



Q2 uvQ'^ - (fc - k'y 



4 



(19) 



with 
/i 

/2 



/a 
3 

/a 
4 



{ap){bp') - {bp){ap) (ug{u) - uvvc 
{ap){bk) — {bp){ak) ug{u) 

{ap'){bk) — {bp'){ak) vvg{v) , 
{ap){bp) — [bp){ap) k"uvvg(v)/2. 
a" {bk) ~ 6" (afc) 1 {pp') u vv g{v) /2 , 
a" [bp] — b" (ap) {pp') uu vvg{v) , 



(v) 



fe = {ap){bp) - {bp){ap') (ap) {v - v) uug{u) , 

= [a" (5p) - fe" (ap)] (py)(ap') uSi;z;g(«) , (20) 

where we have spht the result into different terms to 
facihtate the subsequent discussion. 

In it is understood that the derivatives d/dk" 
act also on the vectors k that are implicit in the func- 
tions /" and in the factors (a/„). Likewise, the ex- 
change of variables indicated in the last line of (ITOl) 
applies also to the functions /j, /" and the factors 
(alu) and (aly). 

One can recognize from the factors g{u) and g{v) in 
((20)) that the hard-scattering graphs with the insertion 
of the chiral-odd operators ([5]) and dH) pick out a twist- 
two distribution amplitude in one of the two pions and 
a twist-three distribution amplitude in the other, as 
anticipated earlier. 



3 Extracting the leading terms 

The factorization formalism is based on an expansion 
in the small parameter A/Q, where A stands for non- 
perturbative momentum scales. In this section we will 
extract the leading terms in this expansion. 

In the following we will assume that the Gegenbauer 
series for g{u) in (1121) converges in the interval u G 
[0, 1], so that 0(it) in (fTT|) vanishes linearly at the end- 
points. The possibility that this may not hold for low 
or moderate factorization scales /i has been discussed 
in a number of papers, see for instance [19l [20l [211 [Ml 
[53]. However, the anomalous dimensions 7^ in ([T5|) 
are positive and increase for n > 0, and evolution to 
high scales will eventually ensure the convergence of 
(fT^ irrespective of the starting conditions. Since we 
are interested in the large-Q'^ behavior, the assumption 
that g{u) is finite at the end-points u = and u = 1 is 
therefore justified. 



Due to the denominators of quark and gluon propa- 
gators, the integrals over u and v in ()16|) can be diver- 
gent when k and k' are zero. From (IT^ and (j^O)) we 
see that these divergences are at most logarithmic in 
both u and v. For the moment we will keep the trans- 
verse momenta k and k' fixed, and regard them as of 
order A Q for the purpose of power counting. We 
first identify terms in (fTO)) that after integration over 
u and V vanish like a power of A/Q (possibly times a 
power of InQ/A). We neglect these terms since other 
contributions will turn out to be finite or to grow like 
a power of InQ/A in the large-Q^ limit. 

To simplify expressions, we use that 



(fk(fk'k°'s{k^,k''^,kk')=0 (21) 



and 



(fk(fk' k"kf^ s{k^,k'^,kk') 

= ^9^ J d^kd^k' k^s{e,k'^kk') (22) 



because of rotational invariance in the transverse 
plane, where s is a scalar function and 



Q/3 

9t 



pp' 



(23) 



Relations analogous to ([2T|) and (|22l) hold with one or 
both of fc", fc^ replaced by /c'", fc''^. 

We now discuss the different terms of ([T^ in turn. 
The reader not interested in the intermediate steps of 
the argument may skip forward to eq. (|55|) . Let us start 
with the contribution involving If the derivative 
d/dk acts on the factors {bk) and (afc) in /f, the result 
is proportional to aag^^bp — bag^^a/B and hence van- 
ishes. If the derivatives act on a factor (afc) in {alu), 
one is left with at least two powers of fc or k' in the 
numerator (a single power giving zero after angular in- 
tegration), which are multiplied by a term proportional 
to 

^ . (24) 

After integration over u and u, this term behaves like 
InQ/A times an even power of A/Q and can hence be 
neglected as well. The terms where the derivative d/dk 
acts on the propagator denominators are proportional 
to 



{ak){bk) ^ {ak){bk) 



aK)[OK) - {ak'){bk) ; , 1 
Mf;g2 - (fc - fc')2 ^ J 

{ak'){bk) - {bk'){ak) 
' uvQ^ - (fc - fc')^ ' 



(25) 
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which vanishes after angular integration. The contri- Furthermore, the contribution from /2 is power sup- 



bution from can hence be neglected altogether. 

We proceed with the contributions from and /" . 
When d/dk acts on a factor (afc) in (alu), we obtain 

{pp') {aa9T'^a/:i){bp) - {ua g^'^b/:s){ap) 

= {ap) [{ap)ibp') - (6p)(ap')] (26) 

multiplied by an expression that, due to the factors u 
and V in the numerator, gives a finite integral over u 
and V even if A: = fc' = 0. If, however, the derivative 
acts on the propagator denominators, we obtain a term 
proportional to 

(ak) 



uQ^ - k^ uvQ'^ -{k- k'f 

[ak) - {ak') 



uvi 



[k - k'Y_ 



fc2 



(6p)-{ao6}. (27) 



At least one more power of (fca) from (aZ„) is required 
to get a nonvanishing term after angular integration. 
The integrals over u and v are only logarithmically 
divergent, so that this contribution is suppressed by 
an even power of A/Q and can again be neglected. 

Let us now discuss the term with f^^. The contri- 
bution from the derivative d/dk acting on fc" needs to 
be retained, whereas contributions with the derivative 
acting on a factor (afc) in {alu) can be neglected: they 
have at least two powers of fc in the numerator, which 
are multiplied by an expression that gives only a log- 
arithm InQ/A after integration over u and v. When 
the derivative acts on the propagator denominators, 
we get a term proportional to 



fc2 



UQ2 

fc(fc 



fc2 

fc') 



UWQ2 _ (yt _ fc/)2 



(28) 



The integral over v of this term gives a logarithm 
InQ/A, whereas the one over u diverges linearly for 
fc = fc' = 0. For finite k and fc' the u-integral thus 
provides a factor l/A^ that cancels the factor A^ from 
the transverse momenta in the numerator. Note, how- 
ever, that the expression in is multiplied by n — 1 
powers of (a^„) = aP{u — ^u) + ak. Only the contri- 
butions from {aP)u need to be retained, since a factor 
u turns the linearly divergent w-integral of (I28|) into a 
logarithmically divergent one, whereas factors of (afc) 
directly provide further powers of (A/Q)^. 

After performing the derivatives d/dk and d/dk' in 
P9p . we can omit all terms (afc) in {aly) and (afc') in 
(aZ„), since they give rise to power suppressed terms. 



pressed and can be neglected. 

Putting everything together we have 



du dv Sk Sk' S(u, fc^) S(w, fc'^) /(u, v] fc, fc') 

{ap){bp') - {bp){ap')\ (aP)--' 

X j dudv d^kd^k' T.{u,k^)T.{v,k''^) 
1 1 f uQ^ 



X 



uvQ"^ - (fc - fc')2 \uQ'^ " fc2 
ug{u) + (u — u) vvg{v) {u — S,u)"~^ 



u"vvg{v) 



k{k~k') 



uQ^ - fc2 uvQ^ - (fc - fc')2 
— (1 + ^) uuvvg{v) {n — l){u — S.u)"'^'^ 

+ (1 + £,){v — v) uug{u) 

^ J2'^u~^u)^-^{v + ^v)"-^-' 
i=i 

— (1 — ^^) uuvvg{v) 

X ^(j - l){u ^ ^uy-^v + ^v)"-^-^ 



+ |ii O w, fc O fc', ^ — > 



In 



91 
A2 



(29) 



Before proceeding let us mention that we checked the 
gauge independence of our result for a general covari- 
ant gauge. Using the same methods as those leading 
to (1^^ , we find that the gauge dependent part of D),^ 
gives only contributions suppressed by an even power 
of A/Q. 

Let us now rewrite (j29l) in a form that allows us to 
identify those terms that give logarithms in Q/A. For 
the term proportional to vvg{v) in the fifth line of (|29p 
we can write 



{u -u){u- ^m)"~i 
= l-2M(u-^t2)"-i 



\Tl-l 



1 - (li - (uY 

= 1 - 2u{u - ^u)"-i -{i+^)uj2{u- ^^y~\ 



(30) 



where in the last step we have used the geometric se- 
ries. Similarly, the terms proportional to ug{u) in (|29|) 
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can be rewritten as 

.1 



n-l 



J = l 
l2 



n-l 



n-l 



- (1 + s - e^)^-' [i - + es)"-^- 

n-l 

- 2(1 + c) ^ (u - (w + iv^-^-^ 



)2 - fc2 



n-l 



n-l-j 



n-l 



1 = 1 



- 2(1 + " + (31) 

In the term proportional to we only need to keep 
the factor u"~^, since with one or more factors of 
we get only a logarithmically divergent integral over 
u and V multiplied by k'^, which is power suppressed. 
Finally, we observe that for those terms in the large 
braces of (EHl) that contain a factor uv, we have 



uvQ^ 



- (A: - k'f uQ^ - P 



l + (32) 



Using the definition ^ of the form factors we then 
obtain 



n-l 



(20 B^niiQ ) =47ras-^^ ^ 



X I dudv d^kd^k' Y.{u,k^)Y.{v,k'^) 



- (1 + vg{v) [n - l)u{u - ^u)""^ 

n-l n— 1— j 

-{l-e)ug{u)Y,{u-iuy-^ ^ {v + ivf-^ 

3=1 1=1 
n-l 

- 2(1 + ug{u)Y,W - ^uy-^ (y + 

-ii-e)v9{v) 

n-l 

X ^(j - l)u(u-^uy-2(w + ^t;)"-l-J 
J = l 

12^ 



= 47ra, 



X S(u,fc'2) 



A2 

^ , dudv d^kd^k' I]{u,k'^) 
2Q2 



{k-k'f 



5(")(1 



vvg{v) + u^^ g{u) 



/c2 



f u"'vvg{v) 
k{k- k') 



iQ^ - fc2 ™g2 - (fc - k'y 



Anas— ^^^"^f^ I dudvugiu) 



1^ 

\n-l 



[n — l)v{v + ^ti) 

n-l 



n-2 



Y^iv + ivy-' 



n-l f 

(1 + Y^u~iur-'-H2{v + £,vy-^ 

.7=1 



+ (1-0 



(j - l)v{v + Cw)-'"^ 

+E(«+^^)'"' 

J=l 1 ^ 



A2 



(33) 



X 



uvQ'^ - (fc - fc')^ 



ug{u) + vvg{v) 



uQ^ - fc2 



?2g2 - /c2 
fc2 



u''^vvg{v) 



fc2 



_uQ2 _ fc2 y^g2 _ (fc _ j^iy 
2vg{v){u-e,uY'-^ 

n-l 

(i+0t'5(f)E("-^^)'"' 
j=i 



where in the last step we have changed the summation 
index j — >■ n — j in the double sum. For the terms 
where the quark and gluon propagators have canceled, 
we performed the integrations over k and k' using the 
normalization condition (|15p for S. 

From p3l) we read off an important result: 

1. The ^ dependent terms of the matrix element ^ 
and thus the form factors i3^„j with i > 2 behave 
like l/Q"* at large Q, up to logarithmic corrections 
from the dependence of as, /Itt and g{u) on the 
renormalization or factorization scale, which one 



7 



should take proportional to . 

These forra factors can be calculated in standard 
coUinear factorization, and the regulating func- 
tions fc^) fc'^) we used in the intermedi- 
ate steps of our calculation have completely dis- 
appeared. The reason for this can be traced back 
to ()19|) . where the only ^ dependence comes from 
the factors {alu) and {alv) and is accompanied by 
factors u 01 V according to p^ . These factors 
suppress the end-point regions and turn out to 
make the u and v integrals finite in the coUinear 
hmit fc = fc' 0. 

2. The form factors Sy„Q involve logarithmically di- 
vergent integrals over u and v in the coUinear limit 
and thus give rise to logarithms of Q/ A if we reg- 
ularize these divergences. 

In the following subsections we shall discuss the two 
cases in turn. 

Before doing so, let us comment on the behavior of 
our result (I33p in the limit of vanishing pion mass. The 
parameter /x^r, which originates from the pion projec- 
tion operator ([7]), is proportional to the chiral conden- 
sate and remains finite in the chiral limit. According 
to ([33|) the form factors therefore vahish like 
in that limit, which is simply due to the factor I/tott 
multiplying them in their definition ([2|) . The pion ma- 
trix element in ^ remains finite in the chiral limit. 
Note finally that when calculating the hard scattering 
we have neglected the quark masses, which are small 
not only compared with Q but also compared with the 
typical values of transverse quark momenta, which we 
have retained in the denominators of propagators to 
avoid divergent integrals. 

3.1 The form factors with i>2 

From p3p one can readily extract the expressions for 
the form factors iJ^^ni with i > 2. The integrals over v 
are elementary, as well as those over u if g{u) is explic- 
itly given as a Gegenbauer series (1121) . For general n 
and k the expressions become rather lengthy, but they 
remain short for the term fc — n — 1 with the maximal 
power of ^. We obtain 



X <; nv ~ -{n + l)v"-^ + 2w"-2 

n-l 



du dv ug{u) 



Ana 



du ug{u) 



n{n — 1) 



2u 



n-2 



E 

i=2 



i-u) 



^ (34) 



where n > 3 must be odd. For n = 3 this gives 



S;^32 = 47ra, 



Cf Gf^m^fi^ . , , , ^ 
gi I duug{u)[-+u 



2 + E' 

n=2 



(35) 



These expressions hold up to power corrections in 
A^/Q^ aj^fj leading order in as- 

3.2 The form factors S^^q 

The form factors Bj-no correspond to the ^- 
independent part of (l33t . Let us first take a closer 
look at terms that have a factor or k{k — k') in the 
numerator. By explicit integration we find that the 
integrals / du dv of 



uvQ"^ - (fc - k'Y uQ^ - /c2 ' 
1 fc2 



uvQ^ - (fc - k')^ uQ^ - k^ uQ^ - fc2 ' 

1 k{k-k') u 



IQ^ ~{k~ k' f uvQ^ - (fc - k')^ uQ^ - 



(36) 



are finite for fc = fc' = 0, as well as the corresponding 
integrals with extra factors of u and v in the numera- 
tor. We thus have 

S?„o(Q') = 8na^^ f^J dudv d^kd^k' 



X S(u,fc2)S(w,fc'2) 



/2n 9{u){l-u^) 



0{l) 



(37) 



where the boldface symbols indicate that we are now 
using a Euclidean scalar product in transverse momen- 
tum space, i.e. {k — k'Y = — (fe — fc')2. Remarkably, the 
r.h.s. of (|37p is independent of n, i.e. the contribution 
enhanced by powers of InQ/A is the same for all n. 
The contribution indicated as 0(1) does not develop 
logarithms of Q/A and depends on n, as is obvious 
from (|33l) . 

To proceed, we replace (fc — fc')2 in (|37|) by a con- 
stant A2 , which thus plays the role of a typical squared 
transverse momentum in the gluon propagator. With 
the normalization condition (jisp for E this replace- 
ment gives 



S(«,fc2)S(z;,fc'2) 



d^kd^k' 



1 



+ (fe-fc')VO^ t25 + A2/Q2 



(38) 
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Clearly, this is an oversimplification since in general 
the average value of {k — k')^ in the integral will de- 
pend on u and v and cannot be described by a single 
constant A^. However, we consider (|38|) as sufficient 
for our purpose, bearing also in mind that even the 
description of the transverse-momentum dependence 
by a single function fc^) is a simplified ansatz, as 
discussed after eq. ([T5|). 

After the replacement (p8|) we can perform the v 
integration in ([BT)) and get 



In stark contrast to the case of in p4[) and 

([35|) . the result (|42|) depends very strongly on the end- 
point behavior of the twist-two pion distribution am- 
plitude 4'{u), or in other words on the higher Gegen- 
bauer coefficients a„ in the expansion (1121) . One can 
expect that Sudakov effects will weaken this depen- 
dence by suppressing the end-points in u, but to in- 
vestigate this is beyond the scope of the present work. 
One should, however, be wary to take the strong end- 
point dependence in (|42l) at face value. 



Cf Gf^m^lJ'^ 



du g{u) 







x(l-s2)lln^^^ + 0(l)^. (39) 



To make the logarithms of Q/A explicit we use that 



^ 1 uQ'^ + K^ 

au — m -T^ 

u A2 



(40) 



and 



du r{u) In 



72Q^ + A^ 
A2 



/•I _ r /_ a^n" 



2 /•! 



ln|^y dur{u) + 0{l) 



(41) 



if r{u) is finite at w = 0. We note that the term of 0(1) 
in (j40|) is equal to tt^/G « 3.3/2, so that one should 
only use our approximation for ln^((52/A2) ^ 3.3. Our 
final result then reads 



5(l)ln^^ -21n . , 
^ ^ A2 A2 



X du 
Jo 



ff(«)-5(l) 



u-1 



+ ugiu) 



Oil) 



24:TTa 



Cf flm^n^ 



)2 



X jln^ (1 + + 15a4 + 28ae + ■■■) 



-In-^ (l + 31a2-l-106a4-l-233.4a6 
A^ 



(42) 



where [g{u) — 5(1)] /(w — 1) is finite at u = 1. 



4 Numerical illustration 

In this section we give some numerical illustrations of 
our results. This is to obtain a basic feeling for the 
order of magnitude and the behavior of our expres- 
sions ([55)) and (HI]). To provide a baseline, we also plot 
the electromagnetic pion form factor, calculated in the 
same approximation as (l35t . i.e. in collinear factoriza- 
tion at leading order in ag-. 



IStTQ!, 



CfJI 
CfJI 



du g(u) 

n=2 



(43) 



At experimentally relevant values of the result (|43l) 
receives important corrections from higher orders in o;,, 
and from various types of power corrections [SJ [HI [Zl HI 
ini [24] . It is natural to expect the same of our result for 
i3p32, and even more so for -B^„o, where the strictly 
collinear framework is not applicable. 

In the following we use the one-loop expression for 
as with riF = 4 active quark flavors and ^qqf, — 
181 MeV. This gives asimr) = 0.33 in agreement 
with extractions of the strong coupling form r de- 
cays [15]. For the quark masses we take the value 
{rriu + md)/2 = 3.79 MeV at the scale fio = 2 GeV 
}12| . which according to ([S]) results in /i^r = 2.57 GeV 
at the same scale. To illustrate the dependence on 
the twist-two distribution amplitude, we take either its 
asymptotic form 0(u) — 6uu or a form with 02 = 0.2 
at /zo = 2 GeV and all other Gegenbauer coefficients 
set to zero. The value of a2 just quoted is close to 
what has been obtained in two recent lattice calcula- 
tions [26] [27]. The one- loop scale dependence of /i^ 
and a„ is given in (jl3l) . in particular one finds that 
/i^(/i) behaves like as(/i)~° '*^ for = 4. 

In Fig. [5] we show our result ([55)) for B^^2 along 
with F-^. We have taken /i^ = fQj. i^j^g renormal- 
ization and factorization scales. For a baseline esti- 
mate this is a natural choice, and we will not explore 
here the more sophisticated options discussed in the 
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literature 124] . We see in the fig ure that 
over an order of magnitude smaller than F-j^ already 
at Q2 = 5GeV2. Of course, the difference between 
these form factors increases with because of their 
different power behavior. We note that both B^^^iQ^) 
and Ft^IQ"^) decrease slightly faster than their nominal 
powers 1/Q^ and l/Q^- This is due to the running of 
CKs, which in the case of i?^32 is more important than 
the increase of /Itt with the factorization scale. We fi- 
nally observe that the dependence on the Gegenbauer 
coefficient 02 is weaker for Bt32 than for _FV, which 
is readily understood from the respective expressions 
dSlD and (gSl). 



0.01 



0.001 







a2(Ho) = 






= 0.2 




10 


100 






[GeV^] 












82(^0) = 0.0 






= 0.2 









10 100 

[GeV^] 



Figure 2: The form factors B'ij.^^ ^^'^ -^V in collinear 
factorization, as given in ((35|) and (|43|) . The factoriza- 
tion and renormalization scales are set to fi = Q. The 
solid (dashed) curve is for a2 — (0.2) at /io = 2 GeV, 
with all other Gegenbauer coefficients set to zero. 



ization and factorization scales may be to take the ge- 
ometric mean /i^ = AQ, which we take as a default 
in the following. In the first panel of Fig. [3] we com- 
pare the results obtained with this choice and with the 
naive choice fi^ = Q^. The differences are noticeable 
but not as large as the ones we discuss next. 

In the second panel of Fig. [3] we compare the form 
factor calculated with three different values of the ef- 
fective parameter A, where the central value A = 
500 MeV corresponds to an estimate based on a model 
of the pion wave function 7J, as discussed in the ap- 
pendix. 

In the third panel of the figure we investigate the 
sensitivity of our result to the twist-two pion distri- 
bution amplitude. The difference between the three 
example choices for the lowest two Gegenbauer coeffi- 
cients are quite small at high but very noticeable 
as decreases. We note that the two curves with 
02(a*o) = 0.2 have a zero crossing, which occurs at 
Q2 ^ 7.8 GeV^ for a^ifio) = and at = 12.2 GeV^ 
for a4(/io) = 0.02. This behavior can be understood 
from (|42p . Compared with the term proportional to 
In^ Q/A, the contribution linear in In Q/ A has a global 
minus sign and larger numerical coefficients multiply- 
ing the a„. If In Q/ A is not large enough, the linear 
term can therefore dominate and give a negative re- 
sult for positive a„. As we discussed after ((42|) . the 
strong enhancement of contributions from higher a„ is 
to taken with great caution, and we therefore do not 
regard the occurrence of a zero crossing for -B^„q as a 
reliable prediction. 

We note that all curves in Fig.[3]fall less steeply than 
a pure power law 1/Q*. This is to be expected since 
the enhancement by the squared logarithm of Q^/A^ is 
stronger than the decrease from the scale dependence 

of as Mtt (m) ^ (m) "'^^ ■ 

Let us finally compare the different form factors 
for our default choices /z^ = AQ with A — 500 MeV 
and a„ = 0. The ratio B^no/ ^T32 varies between 
35 and 240 for between 10 and 1000 GeV^ At 

= 10 GeV^ we find that S^^q is about two thirds 
of Ftt- It is amusing that we obtain i?^„Q = 0.038 at 

= 2.5 GeV'^, which is within a factor of a few from 
the results obtained for iJ^^p ^T20 the lattice 
calculation This coincidence must, however, not 
be over-interpreted, given the uncertainties we have 
just discussed and given that we have not evaluated 
the 0(1) contribution in (|42|) . which is different for 
different n in i3y,jQ. 



Let us now take a look at our result ([^^ for B^^q. 
Since the loop integral in (j37p receives contributions 
from gluon virtualities ranging all the way from order 
to order A^, an adequate choice for the renormal- 



5 Summary 

We have studied the tensor form factors of the pion 
at large squared momentum transfer Q^. The matrix 
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Figure 3: The result (gS]) for the form factor i?r„o- 
Unless specified in the figure keys, we set the renor- 
malization and factorization scale to = AQ with 
A = 500 MeV. As a default we take all Gegenbauer co- 
efficients a„ to be zero; the reference scale for nonzero 
values of a„ is fiQ ~ 2 GeV. 



element of the chiral-odd quark currents with twist 
two are written as the convolution of a hard-scattering 
kernel, the twist-two distribution amplitude for one 
pion and the twist-three distribution amplitudes for 
the other pion. In the twist-three sector we take the 
asymptotic form of the two-particle distribution am- 
plitudes, so that the three-particle distribution ampli- 
tudes do not contribute [TOl [H] . 

For the ^-dependent part of the matrix element ([2|), 
i.e. for the form factors with i > 2, one can take 

the collinear limit of the hard-scattering kernel. The 
result is a representation in standard collinear factor- 
ization, in full analogy with the well-known expres- 
sion (|43| for the electromagnetic pion form factor Ft^. 
The form factors with i > 2 behave like l/Q^ 

up to logarithms from the scale dependence of as and 
/i^ — m'^l{mu + md). Numerically, we find that -6^32 
is more than a factor 10 smaller than F,^ already at 
q2 = 5GeV2. 

For the form factors i?^„Q the collinear limit cannot 
be taken, because the hard-scattering formula then de- 
velops logarithmic divergences in the integrations over 
the longitudinal momentum fraction of the quark in 
both the incoming and outgoing pion. We have used 
a simple regularization of the collinear divergences, 
which involves an effective parameter A representing 
the typical transverse momentum in the gluon propa- 
gator of the graphs in Fig. [T] The momentum fraction 
integrals then give enhancement factors In^Q/A and 
InQ/A that modify the l/Q^ power behavior of B^^q . 
This is reminiscent of the analysis in [28^, where the 
power behavior of the proton Pauli form factor 
F2{Q^) was found to be modified by a squared loga- 
rithm In^ Q /A related with end-point divergences in a 
purely collinear calculation. 

We have evaluated the logarithmically enhanced 
terms for -Bj^rio('3^) ^^'^ ^'^^^ they are indepen- 
dent of the moment index n. These terms depend very 
strongly on the end-point behavior of the twist-two 
distribution amplitude or equivalently on the 

Gegenbauer coefficients a„ with high n. We expect 
this dependence to be decreased by Sudakov effects, 
which suppress the end-points at sufficiently large Q^. 
Numerically, we find that for > 10 GcV^ our ap- 
proximation of i?^„o considerably larger than i?^„j 
with i > 2, which is a direct consequence of the en- 
hancement factor In^ Q/A. 

In the present work we have deduced the basic be- 
havior of the form factors i?^„j at large Q^. An evalua- 
tion that could claim to be quantitatively valid at mod- 
erately large would need to use a formalism with a 
more realistic treatment of the end-point regions in the 
momentum fractions. Obvious candidates for this are 
the modified hard-scattering formalism [6l [71 [HI [22] or 
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approaches based on QCD sum rules [3 HI IS]- 
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A A simple estimate of A 

In order to get some feeling for the typical size of the 
effective parameter A, let us take a closer look at the 
replacement of (fc — fc')^ by A^ in ([55]) . To this end 
we assume that fc^) is independent of u, so that 
we can still perform the integrations over v and u as 
in dMl) to (1121). The logarithms [ln(QVA2)]^ with 
p = 1 , 2 in (H^ should then be replaced by 



j SkSk' s(fc2)s(fc'' 



In- 



{k-k'Y 

Let us for simplicity assume a Gaussian form 



E(fc2) 



271(7' 



■ exp 



'2a2 



(44) 



(45) 



where is the average squared transverse momentum 
in the pion wave function. In a study of using the 
modified hard-scattering picture of Li and Sterman, 
this parameter has been estimated as, a k, 350 MeV in 
conjunction with the twist-two distribution amplitude 

4>{u) = &UU 

With (145]) one can readily perform the integrals (|44)) 
after a change of variables from k and fc' to fc + k' and 
k ~ k' . The result is 



1 



(27r)2a4 
= In 
1 



SkcPk' exp 



+ k' 



2a2 



In 



{k~k'Y 



(27r)2a4 
In 



4e-'ro-2 ' 

(fkd?k' exp 



k^ + k 



'21 



2a2 



In 



[k - k'Y 



4e-Tcr2 



(46) 



where 7 = ^ dxe ^ Inx is Euler's constant. The 
term tt^/G can be neglected in our approximation, so 



that we can consistently identify the first and the sec- 
ond expression in (|46l) with In(QVA^) and ln{Q^/A^), 
respectively. We thus find that with the transverse- 
momentum dependence (|45[) of the pion wave function 
we have A — 2e^''^^a « 1.5cr, which according to the 
above estimate for a corresponds to A « 525 MeV. 
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